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Hastings established exponential decay of correlations for ground states of gapped quantum many- 
body systems. A ground state of a (geometrically) local Hamiltonian with spectral gap e has 
correlation length ^ upper bounded as ^ = 0(l/e). In general this bound cannot be improved. 

Here we study the scaling of the correlation length as a function of the spectral gap in frustration- 
free local Hamiltonians, and we prove a tight bound ^ = 0{l/y/e) in this setting. This highlights a 
fundamental difference between frustration-free and frustrated systems near criticality. The result is 
obtained using an improved version of the combinatorial proof of correlation decay due to Aharonov, 

Arad, Vazirani, and Landau. 


Exponential decay of correlations is a basic feature of 
the ground space in gapped quantum many-body sys¬ 
tems. The setting is as follows. We consider a geomet¬ 
rically local Hamiltonian H which acts on particles of 
constant dimension s; i.e., the Hilbert space is (C®)®", 
where n is the total number of particles. The particles are 
located at the sites of a finite lattice (of some arbitrary 
dimension). We write the Hamiltonian as 

i 

where distinct terms Hi,Hj are supported on distinct 
subsets of particles. Here the support of a term Hi is the 
set of particles on which it acts nontrivially. We assume 
that H has constant range r with respect to the usual dis¬ 
tance function d on the lattice, the shortest path metric. 
This means that the diameter of the support of each term 
Hi is upper bounded by r (e.g., r = 2 for nearest-neighbor 
interactions). Without loss of generality we assume that 
the smallest eigenvalue of each term Hi is equal to zero, 
and that \\Hi\\ < 1. 

If H has a unique ground state |'0) and spectral gap e, 
connected correlation functions decay exponentially as a 
function of distance In particular, 

( 1 ) 

holds with 

e = 0(l/e), (2) 

where d{A, B) denotes the distance between the supports 
of two (arbitrary) local observables A, B, and C is a pos¬ 
itive constant which depends on r and the lattice. In 
the transverse field Ising chain the scaling ^ = 0(l/e) is 
achieved Si , which shows that the upper bound on f 
in Eq. ([2|) cannot be improved. 

In gapped systems with (exactly) degenerate ground 
states, a modification of Eq. dD, 

MAB\^) - {^\AGBM < C\\A\\\\B\\e-‘^^^’^'>/^, (3) 


holds with ^ = 0(l/e) for any ground state {ip) [g, where 
G is the projector onto the ground space. An overview 
of these results and the proof techniques used to obtain 
them is given in Ref. Q. 

Here we specialize to frustration-free geometrically lo¬ 
cal Hamiltonians. Erustration-freeness means that any 
ground state of H is also in the ground space of each 
term Hi. Since we assume that Hi has smallest eigen¬ 
value zero, this means that any ground state |'0) of H 
satisfies Hilfj) = 0 for all i. The ground energy of H 
is therefore zero and in general the ground space may 
be degenerate. The spectral gap e of H is defined to be 
its smallest nonzero eigenvalue. Henceforth we assume 
(without loss of generality Q) that each term Hi in the 
Hamiltonian is a projector, i.e., Hf = Hi. 

Erustration-free Hamiltonians are widely studied in 
physics and quantum computation. 

In quantum complexity theory there is a powerful anal¬ 
ogy between local Hamiltonians and constraint satisfac¬ 
tion problems Q. If we view the terms Hi as constraints, 
then computing the ground energy is a quantum con¬ 
straint satisfaction problem. In this language frustration- 
freeness means that all constraints can be satisfied simul¬ 
taneously. (Classical) satisfiability is a computational 
problem in which one is given a set of constraints and 
asked to determine if they are simultaneously satisfiable. 
Its quantum analogue, defined by Bravyi [lOj (see also 
Refs. [nH3), is a computational problem in which one 
is given a local Hamiltonian and asked to determine if it 
is frustration-free. 

The study of quantum states where the entanglement 
structure is encoded in a tensor network has led to a 
paradigm shift in quantum many-body physics (see, e.g., 
Ref. [hi). There is a close connection between such ten¬ 
sor network states and frustration-free Hamiltonians: for 
matrix product states or projected entangled pair states 
one can construct a frustration-free parent Hamiltonian 
as a sum of local projectors, where each projector an¬ 
nihilates the local reduced density matrix of the state 
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Gapped frustration-free systems include widely stud¬ 
ied spin chains such as the AKLT model [2l| and the 
spin-1/2 ferromagnetic XXZ chain (with kink bound¬ 
ary conditions) |22 h^ . A prevalent strategy to study 
topological phases in two and higher spatial dimensions 
is to construct exactly solvable models, such as the toric 
code [ 2 ^ and, more generally, the quantum double 251 
and string net [i^ models. Almost all such models are 
gapped and frustration-free (commuting, even). Gapless 
frustration-free systems include the spin-1/2 ferromag¬ 
netic Heisenberg chain and the Rokhsar-Kivelson quan¬ 
tum dimer model 27|. Two recent papers construct 
gapless frustration-free spin chains in which the half¬ 
chain entanglement entropy diverges in the thermody¬ 
namic limit n — > 00: a spin-1 example based on paren¬ 
thesized expressions (with log(n) divergent [ 2 ^ and a 
higher spin generalization (i/n divergence) [29|. The clas¬ 
sification of spin-1/2 chains in Ref. [s^ provides many 
further examples of gapped and gapless frustration-free 
systems. 

We establish the following tight upper bound on cor¬ 
relation length in frustration-free systems near criticality 
(i.e., in the limit e —>■ 0). 

Theorem 1. Suppose H is a frustration-free geometri¬ 
cally local Hamiltonian with spectral gap e. Then decay 
of correlations holds with 


C = 0(l/v^) 


(4) 


for any ground state {if) of H. 

In this theorem C in the bound ([3|) is an absolute con¬ 
stant, while the constant hidden in the big-0 notation 
depends only on the interaction range r and a parameter 
g defined by 


g = max|{Hj : yf 0}|. (5) 

I 

{g itself only depends on r and the geometry of the lat¬ 
tice.) 

Theorem [T] may be of interest for at least three rea¬ 
sons. First, for matrix product states or projected entan¬ 
gled pair states, it implies an upper bound on the energy 
gap of the parent Hamiltonian in terms of the correlation 
length of the state. Second, while frustration-free mod¬ 
els such as those discussed above seem to be representa¬ 
tive of many gapped phases of matter, our result states 
that gapless frustration-free systems cannot exhibit crit¬ 
ical phenomena with, e.g., f = 0(l/e). This may be rel¬ 
evant to understanding possible scaling limits of critical 
frustration-free systems, an issue which has been raised 
in Refs. [28l - l^ . The third and final reason, explained 
below, is its potential relevance to the area law for the 
entanglement entropy in one-dimensional spin systems. 

There is a folklore argument which relates correlation 
length and entanglement entropy in one-dimensional spin 


systems. Gombined with Theorem [1] this argument sug¬ 
gests an improved area law for ground-state entangle¬ 
ment in frustration-free gapped one-dimensional systems. 
The area law states that the ground state entanglement 
entropy S for a contiguous region of a (possibly frus¬ 
trated) one-dimensional spin system is upper bounded 
by a constant that is independent of the size of the re¬ 
gion (but which depends on the energy gap e) jsij. In 
particular [^. l^[34| 


S = 0{l/e). 


( 6 ) 


This result agrees with the non-rigorous folklore argu¬ 
ment that particles should be almost uncorrelated if their 
distance is beyond a constant multiple of the correlation 
length. One expects S = 0{f) because only 0{f) par¬ 
ticles in a neighborhood of a cut should contribute non- 
negligibly to the entanglement across the cut. Thus, 
is suggested by ([5]) up to a polylogarithmic prefactor. We 
emphasize that it is at least very challenging (if not im¬ 
possible [ 33 ) to make this argument rigorous. Indeed, 
the connection between entanglement entropy and cor¬ 
relation length has only been proved in a weaker sense: 
S = exp[0(^)] [ 3 ^. Nevertheless it is interesting to note 
that the same argument along with Theorem [T] indicates 
that the stronger bound 


S = 0{l/V~e) 


(7) 


might hold in frustration-free systems. This point, and 
further evidence for Eq. ®, is discussed in Ref. [37|. 

We now turn to the proof of Theorem [1] Reference 
[ 3 ^ gives a combinatorial proof of correlation decay for 
frustration-free Hamiltonians. That proof gives an upper 
bound f = 0(l/e). To prove Theorem [1] we modify the 


argument from Ref. 


using Chebyshev poly nomials. 


In the field of Hamiltonian complexity [39|, l40| , Cheby¬ 
shev polynomials have been used to prove area laws for 
the entanglement entropy in the ground states of one¬ 
dimensional gapped systems 32|, [33|, l4l| . It is thus not 
surprising that they are useful in the present context . 

Proof. —The first part of the proof follows Ref. (38| . 
For completeness, we review the necessary material from 
that paper; we indicate below where this proof differs. 
Here we use a slightly different version of the detectability 
lemma . A proof of this version is given in Ref. [i^ . 


Detectability Lemma ([38|, l43|). Let H = ® 

frustration-free local Hamiltonian with ground space pro¬ 
jector G and spectral gap e. Choose some ordering of the 
terms Hi and let p = n,(i - H.), where the product is 
taken with respect to this ordering. Then, 


\\p-g\\<i/./TT7/? 

where g is given by Eq. ®. 


( 8 ) 
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The result dH) holds for any order of the projectors 
{1 — Hi) in the definition of P. Following Ref. 3^, we 
fix a particular order as follows. It will be useful to define 
an “interaction graph” with a vertex for each term Hi and 
an edge between two vertices if the corresponding terms 
do not commute. Note that g is the maximum degree of 
this interaction graph. 

We first partition the projectors (1 — Hi) into a con¬ 
stant number c of layers (sets) such that any two projec¬ 
tors within a given layer commute. This partition can 
be obtained from a proper vertex coloring of the interac¬ 
tion graph using c colors (no two vertices with the same 
color share an edge). Since any graph with maximum 
degree A has such a coloring with A -|- 1 colors we get 
c < g + 1. For example, for nearest-neighbor interactions 
in one dimension, each projector Hi has support on par¬ 
ticles + 1 , and we may take c = 2 (with one layer 
consisting of all projectors with even values of j, and the 
other layer corresponding to odd values of i). After fix¬ 
ing the layers, we then choose some (arbitrary) ordering 
of them; e.g., in one dimension we might take the odd 
layer to be first and the even layer to be second. Finally, 
we take P = L^ - ■ ■ L 2 L 1 , where Lj is the product of all 
projectors (1 — Hi) in layer j. Choosing P in this way 
we have 13811 


('0|A(P1'P)™B|V') = {iplABlil}) 


(9) 


for m < d{A,B)/{{2c— l)(r — 1)). To see why Eq. (0 
holds, we view (plp)™ as consisting of (2c— l)m layers. 
The reader may find it helpful to look at Fig. 4 in Ref. 
[ 3 ^ . Note that (1 — Pi)]?/') = {ip)] i.e., any projector 
(1 — Pi) acts as the identity on a ground state |'(/'). Like¬ 
wise, for any term Pi with support disjoint from that of 
A we have (1 — Pi)A|^) = \ijj). More generally, in the 
expression (^/>|A(pf P)™ we may replace many of the pro¬ 
jectors with the identity; the ones that remain are said 
to be in the causal cone of A. Each layer reduces the 
distance between B and the causal cone of A by at most 
(r — 1). So if m{2c — l)(r — 1) < d{A,B) then every 
projector in the causal cone of A acts trivially on B\ip) 
and Eq. ([9]) follows. 

At this point we depart from the proof given in Ref. 
[ 3 ^ . using ideas from Ref. [sH. Equation dH) directly 
implies that for any degree-m polynomial Qm{x) with 
Qm(l) = 1 and m < d{A, P)/((2c — l)(r — 1)) we have 

(V'|AQ™(ptp)P|V^) = (V’lAPiV’). (10) 


We choose Qm to be a rescaled and shifted Chebyshev 
polynomial defined by 

T ( _I'l 

Qjn{x) = - 7 -where 6 = , ( 11 ) 

and Tm{x) = cos(marccosa:) = cosh(marccoshx) is the 
standard (degree-m) Chebyshev polynomial of the first 


kind. The function Qm is a degree-m polynomial with 

Qmil) = 1 and [13 

\Qm{x)\ < forO<x<l-A (12) 

Elation HU) follows from Eq. m and the facts that 

\Tmix)\ < 1 for |x| < 1 

Tmix) > for x > 1. 

Since G projects onto the -1-1 eigenspace of the positive 
semidefinite operator plp, we have plP — G > 0. Using 
the detectability lemma we get 

liptp - G|| = IIP - Gf < = 1 - <5 

and, therefore, we have the operator inequality 

0<P^P-G< (l-5)-l. (13) 

Now, let G^ = 1 — G. Again, using the fact that G 
projects onto the -1-1 eigenspace of plp and the fact that 
Qm(l) = 1, we have 

QUP^P) - G = G^{Qm{P^P) - G)G^ 

= G^Qm{P^P-G)G^. (14) 

Using Eq. m and then Eq. d we have, for all 
m < d{A, P)/((2c — l)(r — 1)), 

|(V^|AP|V^)-(V^|AGP|V')| 

= MAiQm{P^P)-G)Bm 

< \\A\\\\B\\\\Qm{P^P)-G\\ 

< \\A\\\\B\\\\Qm{P^P-G)\\ 

< 2||A||||P||exp (^-2m ^, (15) 

where in the last inequality we used Eqs. d and 
d. We now choose m to be the largest integer less 
than d{A, P)/((2c — l)(r — 1)). Substituting the bound 
m > d{A,B)/{{2c— l)(r — 1)) — 1 in Eq. (ITSl) and us¬ 
ing the fact that e/{g'^ + e) < 1, we arrive at the desired 
bound dSl) with G = 2e^ and 

□ 

A simple example shows that the upper bound on ^ 
in Theorem [T] cannot be improved. Consider the spin- 
1/2 ferromagnetic XXZ chain with kink boundary con¬ 
ditions [ 22 I 424 I I . The Hamiltonian for the chain of length 
n can be written as a sum of projectors 

H{q) = |</(<z))(0(g)|.,+i, |0(q)) = 

i=l V9" + 1 

( 16 ) 
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where 0 < g < 1 and |0),|1) are spin up and down, 
respectively. The spectral gap of H{q) is given by [24 1 


e= 1- 


——j- cos(7r/n) 


and vanishes as g —>■ 1 and n —>■ oo. 

The total magnetization M = \ ~ ^i) is con¬ 

served. It can be verified by a direct computation that 


l^i) = 


1 -g^ 
1 — 


1/2 




g^-V"|00...0) 


satisfies H{q)\^i) = 0 and is the unique ground state 
in the symmetry sector where M has eigenvalue 1. Let 
A = |(1 — uf) and B = ^(1 — af) for some j > 1, so that 
d{A,B) = j — 1. Then (i/^i |AB|i/)i) = 0 and so 


\{,l;i\AB\i;,) - {^MGB\i;i)\ 

= l(^i|-4|V'i)(V'i|^l^i)l 

= (17) 

where, in the first equality, we used the fact that A and 
B commute with M. For simplicity we now take the limit 
n —>■ oo. If we suppose Eq. holds for some C and 
then Eq. m implies 

(1 - q^fq^d{A,B) < 


Taking logs on both sides and using the fact that ^ does 
not depend on d{A, B) gives the desired lower bound 

^ ^ -21n(l-0(yi)) ^ 

where in the second step we used the fact that e = 1 — 
2/{q -I- g“^) (in the limit n -A oo). 

Remark .—The XX Z chain also seems to nicely illus¬ 
trate the optimality of the bound ® in the detectability 
lemma. The bound states that 1 — ||P — GH = 11(e), 
while we found using numerical diagonalization that for 
the XXZ chain 116)) 1— ||P—G|| is exactly equal to e, for 
all choices of q and n that we tried (we used the aforemen¬ 
tioned two-layer ordering of projectors in the definition 
of P). Presumably this equality holds for all 0 < g < 1 
and n > 2. 
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